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Abstract 

We study vertex corrections to the leptonic electroweak observables in the 
general MSSM at tan j3 < 35. In particular, we address the question of whether 
supersymmetry can be responsible for the observed 2a deviation from the 
Standard Model prediction in the invisible width of the Z. We find that the 
presence of a light (around 100 GeV) chargino and sleptons hinted by the 
5f^-2 measurements makes the agreement with experiment slightly better and 
improves the electroweak fit. 



1 Introduction 

The recent BNL measurements of the muon anomalous magnetic moment have bolstered 
interest in supersymmetric models [0]. These measurements appear to deviate from Stan- 
dard Model (SM) predictions by 2.6cr p[. A conclusive statement can be made only after 
sufficient statistics have been accumulated and the status of the SM theoretical uncer- 
tainties has been determined unambiguously 0. However, should this deviation persist 
and its error shrink, new physics would be required to explain it. Among the candidate 
models for new physics, supersymmetric models seem most promising ||^. 

In addition to this possible deviation, there are a number of other discrepancies of 
similar size between the SM predictions and the experimental values of the electroweak 
observables. In particular, there is a more long-standing Ab "anomaly" which 
manifests in a 2.7a deviation of the combined left-right asymmetry in Z ^ bb decays 
measured at LEP and SLD from the SM prediction. It has been argued that such a 
discrepancy is unlikely to be a result of a statistical deviation (see e.g. the work of 
Chanowitz in Ref. [^]). The possibility of supersymmetric origin of this "anomaly" will be 
pursued in a subsequent paper. In addition, there is a 2a deviation in the invisible width 
of the Z boson which appears as a deviation of the effective number of neutrinos from 
three: 

= 2.9835 ± 0.0083 . (1) 



Implications of these results for various models of new physics have been considered in 
Refs. 0-0 . In particular, it was found that models with R-parity violating interactions , 
two Higgs doublet models at large tan/3 [0, models with large extra dimensions P, and 
models with an extra gauge U{1)b-3l [0 not only fail to mitigate but in fact exacerbate 
the problem by generating radiative corrections of the "wrong" sign. This observation 
has resulted in stringent constraints on such models. 

In this study we analyze the effect of R-conserving supersymmetry on electroweak 
leptonic observables and, in particular, the invisible width of the Z boson. Motivated by 
the supersymmetric explanation of the BNL g-2 "anomaly" , our study is focused on the 
question whether or not the Z invisible width "anomaly" can be explained with the same 
mechanism. In our analysis we perform a global fit to all relevant electroweak leptonic 
observables such as the R-parameters 
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the left-right asymmetries 
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and the forward-backward asymmetries 

AMI) = \AeAi , (4) 

where /i;^ ^ are the ZIl^fiIl,r couplings and I = e, fi, r. The 2a deviation in Riy/e is 
related to the ~ 2a deviation in ahad = 127rrer/j/m|r|. We remark that Ri,/e is not 
measured directly but rather calculated from the Z-line-shape observables. In principle, 
R^/e could be affected by SUSY contributions to r(Z ^ hadrons); for example, a light 
bottom squark may improve agreement with experiment [^. In this work, we concentrate 
exclusively on the leptonic sector. 

We isolate the effect of the vertex corrections which are sensitive to the lepton/chargino- 
neutralino sector of the MSSM. The oblique corrections are parameterized in our fit but 
not used to constrain the model due to their significant model-dependence, e.g. they de- 
pend sensitively on the Higgs sector, squark masses, etc. In addition, we let as(M^) float 
in our fit since its SM value is extracted from Ri. This strategy has been used previously 
and proven useful in placing generic constraints on complicated models of new physics 
[11-0 • We incorporate the electroweak data reported during summer 2000 conferences in 
our numerical analysis. 

We present general formulae for the vertex corrections in terms of the low-energy 
quantities such as the chargino masses and mixings, left and right slepton masses, etc. 
We then impose the condition of the radiative electroweak symmetry breaking and analyze 
the GUT scale MSSM parameters which improve the electroweak fit. However, we stress 



that our conclusions are independent of the assumptions about the high energy structure 
of the theory and can be formulated purely in terms of the low energy quantities. 

An analysis which addresses a somewhat similar question but with an emphasis on the 



effect of the oblique corrections has recently appeared in Ref. ||TT[. We also find a partial 



overlap with earlier work [11^]. Earlier calculations of one-loop vertex corrections in the 



MSSM may be found, e.g., in Ref. [0. 



The paper is organized as follows. In Section 2 we present our SUSY framework. In 
section 3 we calculate the supersymmetric vertex corrections and study the decoupling 
behavior of the SUSY contributions. In section 4 we discuss the fit and our numeri- 
cal results, and in sectio 5 we make concluding remarks. In the Appendix we list our 
conventions and relevant Passarino-Veltman functions. 



2 Supersymmetric Framework 

We will study supersymmetric models with the following superpotential 

W = -H2QiY^'Uj + HiQiY^'Dj + HdX^Ej - /i^i^2 (5) 
and the high energy scale soft breaking potential 

VsB = + {mlf <l>T<l>a - {B^lH,H, + h.c.) + H,hCe% (6) 

1 
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where (p^^^^ denotes all the scalars of the theory which transform under SU(2) as doublets 
(singlets). We generally allow for nonuniversal gaugino and scalar masses. Note that as 
a result of the SU(2) symmetry different isospin components of the doublets have the 
same soft masses at the high energy scale, whereas there is no similar requirement for the 
singlets. At low energies this degeneracy will be broken by the electroweak effects. 

In what follows we use tan/?, moo, Aa, Mi as input parameters and obtain low energy 
quantities via the MSSM renormalization group equations (RGE) given in Ref. ||14||. We 



also assume radiative electroweak symmetry breaking, i.e. that the magnitude of the fx 
parameter is given (at tree level) by 

2 _ m|^^m|^tan2^ 1 2 
" tan2/5-l 2""^- 

The phase of fi {(f) ^) is an input parameter and is RG-invariant. 

At low energies the charged gauginos and higgsinos mix, leading to the following mass 
matrix (we follow the conventions of Ref. [|I^ ; however, we correct their sign error in the 
superpotential) 



M2 \/2Mw sin (3 

y/2Mw cos p fx 



This matrix is diagonahzed by a biunitary transformation 



U* iVV = diag(m^+, m^+) , (8) 



where U and V are unitary matrices. The mass eigenvalues are defined to be non-negative 
and m + > m+. 

Similarly, for neutralinos we have 



(Ml -MzsmOw cos 13 sin sin 

M2 Mz cos 6w cos /3 — cos 6w sin /3 

— sin 6w cos /3 cos 6w cos /3 — 

V sin 6'vi/ sin [3 —Mz cos 6'vi/ sin [3 —fx / 



This symmetric matrix is diagonahzed by a unitary matrix iV, 

A^* Myfi = diag(m^o, m^o, m^o, m^o) , (9) 

where again the eigenvalues are defined to be non- negative and m^o > m^o > ... The 
chargino and neutralino spinors can be split into the left and right components in the 
usual way: 

Concerning the slepton spectrum, the "left" and "right" charged sleptons also mix at 
low energies. However, their mixing is proportional to the lepton masses and is negligible 
unless tan /3 is very large. Neglecting lepton masses, the low energy mass eigenstates are: 



m 
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mf + M| ( -- + sin^ Ow ) cos 2/? 



ml c::^ ml — M| sin^ 9w cos 2/? 



ml ~ + cos 2/3 , (11) 



where m~ and mj are the mass parameters appearing in the low energy analog of Eq. 



3 SUSY Vertex Corrections 

In this paper we will concentrate on tan (3 < 35. It is quite difficult to achieve radiative 
electroweak symmetry breaking for greater tan f5, so such an assumption can be justified. 
At tan P <35 the gauge couplings dominate the lepton Yukawa couplings, so only the 
gaugino parts of the charginos and neutralinos couple to leptons with an appreciable 
strength. In addition, one can neglect the left-right slepton mixing in this regime (as will 
be clear below, each relevant diagram would require two left-right mass insertions, so the 
effect of this mixing is further suppressed). 



We perform our calculations using the two-component spinor technique (see the Ap- 
pendix for the notation and conventions). The result is expressed as a correction Shf^^ 
to the tree level coupling hf^ ^ defined by 
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(12) 



with 



hf^ = I3 ~ Q sin^ 9 
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Neglecting the lepton Yukawa couplings, we have the following SUSY interactionsQ|T5| 

^ -9[{Unxt + U;,^)PLyil + {V*,W' + V*,W')PLeV*] +h.c., 
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Here J3 and Q are the lepton isospin and charge, respectively, and the superscript c stands 
for a charge conjugated spinor. The vertex structures Oij are given by 



-VaV;,--V,2V*, + 6,,sm^9 
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-UnUji - -U*^Uj2 + 5ij sin^ 9w , 
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(15) 



These interactions are to be expressed in terms of the two-component spinors. The im- 
plementation is trivial for all interactions except for which becomes 
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{UnxtR + U2ixtR)eL + eiia2{Vn{xtLr + V2i{xtLr)i>L] + h.c. (16) 



We remark that denotes a Dirac spinor with a positive charge (not to be confused 
with a hermitian conjugated spinor). 

*This corrects an error in Ref. in the expression for the neutrahno couphng to right-handed leptons 
(C77), i.e. N*2 should be N*^. 



The Z-x°-x° coupling can be simplified by taking advantage of the Majorana nature 
of the neutralino. For Majorana spinors ipi and tp2 we have 



Using this identity as well as 0"J^ = —0"jl, we obtain 



(17) 
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(18) 



3.1 Chargino Contributions 

In this subsection we list expressions for Feynman diagrams containing charginos in the 
loop. Since the higgsino coupling to leptons can be neglected at tan [3 < 35, the charginos 
induce corrections to the left-handed couplings only. Below we present our results in 
terms of the corrections to the tree level Z-fi-fL couplings /i/^. 
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(la) : g"" J2 OiV*Vj, [{2 - d)C2A + MfCasJ (M|; m,, m^+, m^p 



(Ic) : -g'^^\Vki\^ (724(^1; m^+,m£;,m^) , 

k ^ 

{Id) ■ [-^- + siYiHw^Y.\^ki? 5i(0; m^+,m,) . 



(19) 



Definitions of the B and C functions can be found in Appendix B. We note that in 
addition to Fig. Id there is another wave function renormalization diagram with the 
loop on the outgoing electron leg. The corresponding correction is the same as for the 
diagram in Fig. Id, so we do not list it separately. The contribution of the wave function 
renormalization diagrams to the total correction comes with a factor of 1/2, so in effect 
the total correction is simply given by a sum of individual contributions in Eq.([TP|). The 
analogous contribution to the (left-handed) neutrino final state is 



5K: 
(2a) 

(26) 

(2c) 



-9^ E 0\fU*^Ua [(2 - d)C2A + MfC-ss] (M|; mg,, m^+, m^+) 
ij 

ij 

-g'^{-l + 2sm^ 9w) E l^^il^ (724(^1; m^+, mg^, mg^) , 
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The resulting total corrections are 



5h' 



Sh'„ 



k ^ 

E 0\fU*,Ua [(2 - d)C'24 + MlC'23] (M|; me,, m^+, m^+) 



+ E ^it^ii^ii "^x+"^x+'^o(^^^' "^^'i' "^x+' "^x+) 



+ (-l + 2sin2^w')^|[/fci|2 (724(^^1; m^+, me,, me J 



(21) 



(22) 



These corrections are finite as they should be. This can be seen from the relations 

E 0'iV*,V^, = E 0[fU*,Ua = -1 + sin^ Ow , 
T.y:lV^k = EUnU^k^S^k (23) 



and the fact div(C'24) — —1/2 div(Si) while Cq and C23 are finite. 
3.2 Neutralino Contributions 

Because of their bino component, neutralinos induce corrections to both the left and right 
couplings of the leptons. Starting with the correction to the right-handed charged lepton 
coupling, we have 

(3a) : -2/ tan^ Ow E 0'-l'N*^Nj, [(2 - d)C'24 + MlC'23] (M|; mg^, m^o, m^o) , 

ij 

(36) : isjn^ 6w^^0^^N*^Nii m^om^p Co{Mz; 171^1^,171^0,171^0) , 



(3c) : -ig"^ tan^ 9w sin^ 61^/ E l^^^il^ C'24(M|; m^o, mg^, m^^) , 
(3d) : —2g^ tan^ sin^ 9w E I -^fci I ^ -Bi(0; m^o, mg^) . 



(24) 



The corrections to the left-handed charged lepton couphng are given by 
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[(2 - d)C24 + MlCsal (M|; mg^, m^o, m^o) , 
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(4c) : (^-^ +sin2^w^^ \Nk2 + tein0wNki\'^ C'24(M|; m^o, mg^ 
(4d) : -|- (^-i + sin' Ow^ Yl \^k2 + tan 9wNki\^ Bi{0; m^o, mgj . 



Finally, the neutrino coupling corrections are 



[(2 - d)C2i + MfCss] (M|; m^, m^o, m^o) , 



(56) : ^-Y.O'-^'iN*^ - tan^^TV* )(Ar.2 - tan^^TV-,) 



X m^om^o C'o(M|; mjj, m^o, m^q) , 
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(5c) : -— ^ |A^fe2 - tan6lw^A^a|^ (724(^1; m^o, mp, m^) , 



The total corrections are given by 



-2g^ tan' Qw 
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These expressions are finite due to the relations 



E^y ^4^.^ = , (fc,/ = l,2) 



(30) 



and the fact that the combination 2(724 + Bi is finite. Note that the diagrams in Figs. 
3a, 4a, 5a are individually finite. The reason is transparent in the weak eigenstates basis: 
only the higgsinos couple to Z, and we retain only the gaugino coupling to the leptons, 
so a mass insertion is necessary on each fermion line to complete the diagram. 

3.3 Decoupling of Heavy Superpartners 

In this subsection we demonstrate explicitly the decoupling of heavy SUSY particles. As 
the SUSY mass scale increases, the gauginos and higgsinos become approximate mass 
eigenstates and V, U can be chosen such that 
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In the expressions for the vertex structures O^^-, the factors in front of the Kronecker delta 
symbols represent the gaugino and higgsino couplings to the Z boson. It is clear that 
corresponds to the gaugino (W^~) Z coupling, while O22 corresponds to that of the 
higgsino {h~). Since O^^ is to be contracted with Vjt^i, the higgsino component drops 



out of all expressions in the decoupling limit, as expected. Denoting by m a heavy scalar 
mass and by M a heavy fermion mass, we can rewrite as 



5h: 



sin' Owi {2 - d)C24 + MI;C23 (0; m, M, M) + M'Co(0; m, M, M) 
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(32) 



Each of the expressions in the curly brackets vanishes, see Appendix B. Note that even 
though the corrections in Eq.|3T| are linear in Mz/mgusy^ the SUSY contributions decouple 
quadratically, as they should. Similarly, for the neutrino final state we have 
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(33) 



Concerning the neutralino contributions, let us first consider ^h!^^. Since the mixing 



between the gauginos and higgsinos vanishes in the decoupling limit, we have 

' Mz \ 
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As a result, the combination O'^^ N*-^N ji vanishes in this limit. Therefore 
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(34) 



. 



Again, the combination in the curly brackets vanishes (see Appendix B) . The same argu- 
ments are valid for the neutralino corrections to the couplings of the left-handed leptons. 



4 Numerical Analysis 



To separate out the effect of vertex corrections, we pursue the strategy of Refs.|l6|-|9[],p^ . 
That is, we utilize only those observables which can be expressed as ratios of the weak 
couplings. The effect of oblique corrections [|17| then either cancels in the ratios or can 



be absorbed into effective sin^^vK- In the fit, we leave siia^ 9w as a free parameter and 
parameterize the vertex corrections as 6h^, 6hi^, and Shi^. In addition, we retain as{Mz) 
as a free parameter since its value is determined from Ri. The fit value of 5(sin^ 9w) is 
not used for constraining the model due to its model dependence. Specifically, (5(sin^ 9]y) 
depends on the Higgs, squark, etc. masses and thus is not particularly useful in our 
general analysis. 

We impose the following (direct search) constraints on the SUSY spectrum p 



me > 99 GeV , 

> 96 GeV , 
rrif > 87 GeV , 

> 43 GeV , 
m^o > 36 GeV , 

> 94 GeV . (35) 

We assume that the lepton parameters are generation-independent since lepton-universality 
breaking corrections are quite constrained (see, for example, the second reference in 0); 
in any case this assumption is not important for our analysis. 

Before we proceed, a few comments are in order. First, note that fi is determined by 
a particular combination of the Higgs mass parameters (i.e. m|^^ — m^^^tan^/? ). Thus 
having fixed fi, the squark masses and the "orthogonal" combination of the Higgs mass 
parameters remain free. This freedom results in the uncertainty in the oblique corrections 
mentioned above. However, to be specific, we will fix them in our numerical analysis still 
assuming the freedom in the oblique corrections. Second, our results are presented in terms 
of the high energy parameters. Since we do not generally assume a particular framework 
or relations among the soft breaking parameters, one might wonder why not interpret our 
results directly in terms of the low energy quantities. However, not every low energy set of 
parameters can result from some high energy boundary conditions, especially consistent 
with radiative electroweak symmetry breaking. To cite just one example, heavy gluinos 
and light squarks at low energies are inconsistent with high energy boundary conditions 



if we are to avoid color breaking minima |l9[. Thus, to be safe, we will generate each low 



energy set of parameters via the RG running. 

Numerically Ru/e is sensitive to the vertex corrections 6hi and much less sensitive to 
the oblique corrections: 

6R^/e = 7.96 5K + 8.50 5he^ - 7.33 + 1.17 5s^ , (36) 

where = siia^Ow- Similarly, for the left-right asymmetries we have 

5Ae = -3.64 5he^ - 4.23 5/ie« - 7.87 5s^ . (37) 

Since the SM prediction for R^/^ is above the measured value whereas that for the lepton 
asymmetries is below the measured values, 6hej^ < is favored by both R^/^ and Ai, ApB- 
As shown below, 5he^ in the MSSM is typically larger than dh^ and dhe^,. To get a feeling 





A'h = 150 


M2 = 170 


M2 = 190 


M2 = 210 


5as X 10^ 
5s^ X 10^ 


11.09 
3.9 ±3.9 
-9.9 ±2.0 


11.30 
3.4 ±3.9 
-10.1 ±2.0 


11.40 
3.1 ±3.9 
-10.2 ±2.0 


11.48 
2.8 ±3.9 
-10.3 ±2.0 



Table 1: The quality of the fit ^is a function of the GUT scale parameter M2 (GeV). 
For comparison, the Standard Model gives = 11.44 (d.o.f.=12-2). The corresponding 
fit values of as and sin^ 9w are also displayed. The other GUT scale parameters are 
= 10 GeV, mg = 85 GeV, Mi = 100 GeV, M3 = 150 GeV, A = 100 GeV, tan/5 = 5 
and the scalar mass parameter (except for sleptons) is set to 100 GeV. 

for the value for preferred by the fit, set 6hi, = Shg^^ = and fit Ri, Ai, Aps with 
three parameters: 6hej^, Sag, and The best-fit values are 

She^ = -0.00165 ± 0.00096 , 
5as = 0.024 ±0.014 , 

^ _0.0002 ± 0.0005. (38) 

R^/e strongly pulls to be negative, resulting in a large correction to Ri which is in 
turn compensated by a large 6as. In addition to a genuine shift in as, our "effective" 
6as parametrizes potential corrections to r(Z — > hadrons) from the squark/Higgs sectors. 
This is, of course, just a "toy" fit. As we will see below, for a viable MSSM model all the 
shifts are much smaller. Qualitatively, however, the picture remains the same - a negative 
She]^ is preferred by the data. 

Consider now the lepton vertex corrections in the MSSM. In Figs.|^-P we display the 
vertex corrections Sh,^, 6hf.j^, and Sh^j^ as functions of M2 and tan/3. Note that Sh,^ and 
are quite sensitive to M2, whereas its effect on is negligible as it arises only via 
the RG running. In most of the parameter space, 6hej^ dominates the other corrections; it 
has the right sign (negative) to mitigate the invisible width "anomaly" , especially for the 
positive sign of the /i-term (which is also preferred by the Qfi — 2 measurement). We find 
that the regions of the parameter space where SUSY contributions improve the agreement 
with the measured values of — 2 and Ru/e are generally compatible, see for instance 



Ref. 20 



In Figs. |T0| and 0, we display the corresponding shifts in R^y/e as functions of M2 and 
tan/? (keeping sin^6'iy fixed). Varying M2 from 135 to 250 GeV corresponds to varying 
the light chargino mass from 95 to 180 GeV. Fig. ^ shows the dependence of Ru/e on the 
GUT left slepton mass parameter mf its range 10-200 GeV translates into the slepton 
mass rngj^ range of 104-225 GeV. 

The error bar for R^/e is 0.008 (see Table 2), so the supersymmetric contributions 
can only be responsible for the shift of about 0.2a. This suppression results partly from 
the cancellation between the neutrino and left-handed electron contributions. Indeed, 
if the chargino is a pure gaugino and the sneutrino and left selectron masses are equal, 
Shu = ~^h(.^ in the Mz — > approximation and the resulting contribution to Ru/e is very 



small (Eq.^) . Alternately, if the chargino is a pure higgsino, the corresponding couplings 
are very much suppressed and the resulting SRyje is negligible. One thus expects the 
largest correction when there is a large splitting between the sneutrino and left selectron 
masses (which is severely bounded by the SU(2) symmetry) and/or when the chargino is 
a gaugino-higgsino mixture {M2I n ~ 1) Q- 

The dependence on other input parameters is significantly weaker. An increase in 
M3 affects the /x-term via the radiative EW symmetry breaking condition, which in turn 
results in heavier charginos and neutralinos. The effect of Mi is not significant due to the 
subdominant role of the neutralino contributions. For the same reason the dependence 
on the masses of the right sleptons is weak. 

For completeness, below we provide representative low energy parameters for our stud- 
ies. The GUT scale parameters mi = 10 GeV, mg = 85 GeV, rriQ = 100 GeV, Mi = 100 
GeV, M2 = 135 GeV, M3 = 200 GeV, A = 100 GeV, 0^ = 0, where mo is the mass 
parameter for the scalars other than sleptons, lead to the following low energy spectrum 

~ (372, 95) GeV , m^o ~ (374, 353, 96, 38) GeV , 
m/^ 75 GeV , mg^ ~ 104 GeV , mg^ ~ 101 GeV , 

and the following mixing matrices 

^ _ . 0.34 0.94 \ ^ /0.19 0.98 



N 



0.94 -0.34; ' V0.98 -0.19 

f 0.11 -0.26 -0.67 0.69 

0.05i -0.08i -0.70i -0.71i 

0.15 0.95 -0.22 0.12 

V-0.98 0.11 -0.15 0.06 / 



(39) 



We now turn to the discussion of the fit. In Table 1 we present our fit results for 
different values of M2. That is, we fix the lepton vertex corrections using our GUT 
parameters in the fit and calculate the corresponding (d.o.f.= 12-2). The parameters 
Ss"^ and (effective) Sag are left as free fit parameters, which means that the Higgs and the 
squark sectors "adjust" themselves so as to give the best fit results. 

To determine if there is any improvement over the SM, we perform a fit for the SM 
under the same circumstances, i.e. lepton vertex corrections set to zero, 63"^ and 6as free 
to account for a variation in the Higgs mass and a^. The Standard Model fit gives 

= 11.44 (d.o.f. = 12 - 2) , 
Sas = 0.0020 ±0.0039 , 

Ss^ = -0.00103 ± 0.00020 . (40) 

If the chargino is light (100 GeV), the MSSM fit gives = 11-09. We see that the SUSY 
vertex corrections indeed improve the fit due to the improvement in R^/e and the lepton 
asymmetries. The quality of the fit quickly approaches that of the SM as the chargino 
mass increases. We note that the best fit value of for both the SM and the MSSM 
significantly deviates from zero because of the SLD asymmetries, which signifies that the 
light Higgs is preferred. 

'^This was also noted in Ref. ImJ 



5 Conclusions 



We have analyzed leptonic electroweak observables in the general MSSM. We find that 
supersymmetry can mitigate some of the discrepancies between the Standard Model pre- 
dictions and the observed values. Namely, it produces vertex corrections of the right sign 
to improve agreement with Ry/e and the leptonic asymmetries. As a result, the elec- 
troweak fit is improved from = 11.44 (SM) to = 11.09 (MSSM). This required a 
light (~ 100 GeV) chargino and relatively light (100-250 GeV) sleptons. 

Although the improvement from statistical point of view is not very significant, it is 
quite encouraging since in the same region of the parameter space the — 2 discrepancy is 
also mitigated. This is to be contrasted with a number of "new physics" models considered 
earlier [^-0, all of which made the electroweak fit worse. The improvement of the fit 
requires light superpartners which can be detected in collider experiments in the near 
future. 
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A Appendix: Notation and Conventions 

We use the following (chiral) representation of the Dirac matrices: 

where = (1, 'd^) and = (1, — "o^). The left and right components of a Dirac spinor 
and the corresponding projectors are defined by 

The charge conjugated spinor is given by 

= , C = -i-i^-i^ , = ^^^7° . (43) 
In terms of the two-component spinors this corresponds to 

Free fermions satisfy the following Dirac equation in the two-component notation 

{k ■ a) = m , 

{k ■ a) i)R = m ipL , (45) 

The corresponding propagators read 

I k ■ a + k ■ a 

m ^i) = i-Hi — , {'4'R Wr) 



f^2 _ yj^2 ' ^2 _ j^2 ' 

m 



{i^L V^L) = {^R ^l) = ^^rzr^2 ■ (46) 



The following identities are useful for calculating Feynman diagrams in terms of the two- 
component spinors: 



2 T 2 — 2 — T 2 



a^ = {2-d) a 



5 



{p ■ a) (p ■ a) = -v" . (47) 



B Feynman Integrals 

Here we make explicit our notation for the scalar and tensor integrals that appear in the 



calculation. The definitions of the integrals are slightly different from those of Ref. 
The hat on the tensor integrals serves as a reminder of these differences. 



B.l Scalar Integrals 

We define the functions Bq and by: 



5o[p^;mi,m2] = 

C'o (P - g)^;'^i,"i2,"^3] = 

The general form of Bq is given by 

-1 



% /i 



4-(i 



1 



{2t^Y (P - m?) [ (A; + p)2 - mi ] ' 
d^k 1 



(4J 



(27r)4 (P - m\) [ (fc + p)2 - mi ] [ (A; + - mi ] ' 



i?o[p';mi,m2: 



(47r)^ 



_ m^ln(m^/;.2)-miln(mi//.^) ^ ^ ^ ^^^^ 
— m2 



where = — 7s + ln47r, and 



2 Vl-S- Vl-2S + A2^ 



(49) 



with 



^ ^ m| + m^ ^ ^ 



mt — ms 



(50) 

The function F{p'^\mi,m2) vanishes in the limit —>■ 0. The general form of the Cq 
function is fairly complex and we refer the reader to Ref. |^T[. The special case relevant 
for our calculations is 



Co [ 0,0,(5^; mi, m2,m3' 



dy 



1 



(47r)2 Jo mi — ml — yQ"^ 



X In 



yjy - l)Q^ + (ml - ml)y + mj 
y{m2 — mf) + mf 



(51) 



B.2 Tensor Integrals 

Definition and general form of Bi: 

1 1 / 771^ Tfi^ \ 

Bi[p'^;mi,m2] = --5o [p^ mi, ] + -^^^ f ^^^^ j mi, mg). (52) 

Note the following useful relations among the i?-functions: 

= Bo[p';m„m,]+B,[p';m,,m,] + B,[p';m„m,], (53) 
= (m^ - mg) -Bo [0; mi, + (m2 - mg) i?o [0;m2,m3] + (mg - m^) i?o [0;m3,mi ]. 

Definition of the C-functions: (Note the difference from the definitions in Ref. |pi|| .) 

C^[p,g;mi,m2,m3] = t J ^^^^^ _ ^ ^^2-^ p + _ ^2] 

= P^,Cii + q^Ci2, (54) 

\ 1 • 4-d f ^'^^ k^ku 
p,g;mi,m2,m3 = TTr^JT^ 2\ \ii , — V2 21 \n , — ^2 2T 

For the purpose of this paper, we will only need to evaluate these functions for = = 
(we neglect final state fermion masses). = {p — qY — ~'^P ' 1 ^^^^ then be the invariant 
mass squared of the initial vector boson. For this parameter choice, the C-functions can 
be expressed in terms of the 5-functions and Cq as: 

Cii = {^o [0;mi,m2] - Bq [Q'^; 1712,1713] - (m^ -m^)C'o} , 
C12 = {^o [0;mi,m3] - Bq [Q^; 1712,1713] - {ml - ml) Co] , 



C 

C24-- 



-Bi [Q2;m2,m3] + {mj - ml) Cu + m^ Cq - 



{d - 2) (724 - ^23 = -fil [ g'; m3, m2 ] - {m\ - ml) Cn + . (55) 

We do not list expressions for 6*21 nor C22 since we do not use them in this paper. 
B.3 Decoupling Limit 

Below we list approximate formulas valid in the decoupling limit p^ /ml ^ 0. Here mg 
and mj denote the scalar and fermion masses, respectively. Omitting the 0{p'^/m'^ y) 



terms, we have 

[{d - 2) C24 - C23] (0, 0,p^; nis, ruf, ruf^ 

C24(p,0,P^;rnf,ms,ms^ R: 
mjCo (0,0,p^; ms,m/,m/) 
Bi (0;m/,ms) 



(Att) 



m. 



1 ( A,-ln-^ 



2(47r)2 
1 



A, - In 



ml 



(An) 



[f{x)+g{x)], 
i('A.-ln^l-5(x) 



where 



4(l-a;)2V /' 
1 , 1 



-- Inx + - , 
2 4(1 



(l-x)(l-3x) + 2x^ Inx 



(56) 



(57) 



for X = mj/rrig. 
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Observable 


Measured Value 


ZFITTER Prediction 


Z lineshape variables 




mz 


Qi 1 S7fi -1- n nnoi HcAT 
y i. io ( It u.uu/i vorev 


1 Tl Til 1 "f" 
lliJJ U-lj 


Tz 


2.4952 ± 0.0023 GeV 


n Tin cdH 


"ha.d 


41.541 ± 0.037 nb 


llTlllSPrl 


Re 


20.804 ±0.050 


20.739 


R/j. 


20.785 ± 0.033 


20 739 


R-j- 


20.764 ±0.045 


20.786 


^FB(e) 


0.0145 ±0.0025 


0.0152 


Afb(/U) 


0.0169 ±0.0013 


0.0152 


^FB(r) 


0.0188 ±0.0017 


0.0152 


Ru/e 


1.9755 ±0.0080 


1 9916 


T polarization at LEP 






0.1498 ±0.0048 


0.1423 


Ar 


0.1439 ±0.0042 


0.1424 


SLD left-rig 


lit asymmetries 




Alr 


0.1514 ±0.0022 


0.1423 


Ae 


0.1544 ±0.0060 


0.1423 


A, 


0.142 ±0.015 


0.1423 


Ar 


0.136 ±0.015 


0.1424 



Table 2: LEP/SLD observables and their Standard Model predictions. The data are from 
Refs. and [^]. The Standard Model predictions were calculated using ZFITTER v. 6. 21 
[|1 with mt = 174.3 GeV, = 300 GeV, and a,(mz) = 0.120 as input. 





T^Z ^Z fT^ad 


Re 


Rfi 


Rt 


^FB(e) 


^fb(/w) 


^FB(r) 


mz 


1.000 -0.008 -0.050 


0.073 


0.001 


0.002 


-0.015 


0.046 


0.034 


Tz 


1.000 -0.284 


-0.006 


0.008 


0.000 


-0.002 


0.002 


-0.003 


'-'had 


1.000 


0.109 


0.137 


0.100 


0.008 


0.001 


0.007 


Re 




1.000 


0.070 


0.044 


-0.356 


0.023 


0.016 


R/i 






1.000 


0.072 


0.005 


0.006 


0.004 


Rr 








1.000 


0.003 


-0.003 


0.010 


^FB(e) 










1.000 


-0.026 


-0.020 


^FB(yU) 












1.000 


0.045 


^FB(r) 














1.000 



Table 3: The correlation of the Z lineshape variables at LEP. The correlation of R^/e with 
Apb{g) is ±0.28, while its correlation with the /i and r observables is negligible. 




(a) (b) (c) (d) 



Figure 1: Chargino-sneutrino corrections to the Zej^j^ vertex (in this and other figures 
the wave function renormahzation diagram with a loop on the upper fermion leg is not 
shown) . 




(a) (b) (c) (d) 



Figure 2: Chargino-selectron corrections to the Zv^vl vertex. 




Figure 4: Neutralino-selectron corrections to the Zei^e^ vertex. 




Figure 5: Neutralino-sneutrino corrections to the Zvi^vj^ vertex. 
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Figure 6: Vertex corrections to the Zff couphngs as a function of the GUT scale 
parameter M2 for 0^ = 0. 1 - 6hy^ 2 - ^/ig^,, 3 - ^/le^- The other GUT scale parameters 
are m^- = 10 GeV, mg = 85 GeV, Mi = 100 GeV, M3 = 200 GeV, A = 100 GeV, 
tan /3 = 3. The other scalar mass parameters are set to 100 GeV and the CP-phases are 
set to zero. 
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Figure 7: Vertex corrections to the Zff couplings as a function of the GUT scale 
parameter M2 for 0^ = vr. 1 - 6hu, 2 - 6hej^, 3 - ^/le^- The other GUT scale parameters 
are m^- = 10 GeV, mg = 85 GeV, Mi = 100 GeV, M3 = 200 GeV, A = 100 GeV, 
tan P = 3. The other scalar mass parameters are set to 100 GeV and the CP-phases are 
set to zero. 
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Figure 8: Vertex corrections to the Z f f couplings as a function of tan (3 for 0^ = 0. 1 
- Shy, 2 - 5hf,j^, 3 - The other GUT scale parameters are m^- = 10 GeV, = 85 

GeV, Ml = 100 GeV, Ms = 135 GeV, M3 = 200 GeV, A = 100 GeV. The other scalar 
mass parameters are set to 100 GeV and the CP-phases are set to zero. 
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Figure 9: Vertex corrections to the Zff couplings as a function of tan (3 for 0^ = vr. 1 
- Shy, 2 - 6hej^, 3 - 6hej^. The other GUT scale parameters are mj = 10 GeV, mg = 85 
GeV, Ml = 100 GeV, M2 = 135 GeV, M3 = 200 GeV, A = 100 GeV. The other scalar 
mass parameters are set to 100 GeV and the CP-phases are set to zero. 
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Figure 10: Shift in R^jf. due to the vertex corrections as a function of M2 (this corresponds 
to the range of the hght chargino mass from 95 to 180 GeV). 1 - 0^ = 0, 2 - 0^ = vr. The 
GUT scale parameters are m,- = 10 GeV, mg = 85 GeV, Mi = 100 GeV, M3 = 200 GeV, 
A = 100 GeV, tan/3 = 3. 
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Figure 11: Shift in R^/^ due to the vertex corrections as a function of tan/5. 1 - 0^ = 0, 
2 - 0^ = TT. The GUT scale parameters are mj = 10 GeV, rrig = 85 GeV, Mi = 100 GeV, 
M2 = 135 GeV, M3 = 200 GeV, A = 100 GeV. 
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Figure 12: Shift in Ry/e due to the vertex corrections as a function of the GUT scale 
slepton mass parameter m^- (this corresponds to the range of the slepton mass mg^ from 
104 to 225 GeV). 1 - 0^ = 0, 2 - 0^ = vr. The other GUT scale parameters are mg = 85 
GeV, Ml = 100 GeV, Ms = 135 GeV, M3 = 200 GeV, A = 100 GeV, tan/3 = 3. 



